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ABSTRACT 

Dispersion  and  impedance  properties  in  various  structural  wave- 
guides, including  non-uniform  rods  and  beams  were  examined.  Here  are 
some  conclusions:  (1)  Even  though  both  a dilatational  and  a shear 
potential  cure  required  to  express  longitudinal  wave  motion  in  3- 

dimensional  cylinders,  the  axial  modal  impedance  equals  pc  , where  p 

z 

is  density  and  c the  modal  phase  velocity.  (2)  At  cut-off  the  axial 
z 

modal  and  radial  modal  impedances  respectively  diverge  and  vanish. 

(3)  Except  for  torsional  waveguides,  whose  wave  motion  is  expressible 
in  terms  of  a single  potential  the  drivepoint  impedance  cam  not  be 
calculated  as  an  eigenmode  series,  as  these  do  not  form  a complete 
orthogonal  set.  (4)  In  non-uniform  longitudinal  waveguides,  increasing 
density  and  Young's  modulus  give  rise  to  attenuation  by  backscattering. 
Whether  these  parameters  vary  exponentially  or  as  a power  of  the  axial 
coordinate,  their  contribution  to  the  attenuation  is  precisely  half  the 
spreading  loss  due  to  a similarly  varying  cross  section.  (5)  Cut-off 
behavior  of  the  fundamental  longitudinal  mode  exists  for  exponentially 
varying  parameters  if  Young's  modulus  and  density  are  identical 
functions  of  the  axial  coordinate;  when  these  parameters  vary  as  powers 
of  the  coordinate,  certain  powers  lead  to  non-propagating  deflections 
or  exponentially  decaying  waves;  wedges  or  cones  undergoing  flexural 
vibrations  display  a cut-off  frequency  when  supported  on  continuous 
springs. 
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DISPERSION  AND  CUT-OFF  PHENOMENA  IN  RODS  AND  BEAMS 


tntroduction  and  Principal  Results  of  this  Studv 


1 


f. 

1 1 


This  study  was  proposed  to  ONR  Structural  Mechanics  Program  as  a 
result  of  in-house  research  which  led  to  the  conclusion  that  uniform  struc- 
tural waveguides  of  small  transverse  dimensions  (1)  display  a varnishing 
impedance  at  every  modal  cut-off  frequency;  (2)  that  each  cut-off  frequency 
coincides  with  the  natural  frequency  of  the  corresponding  standing-wave 
node  of  the  cross  section,  the  modal  configuration  being  independent  of 
the  z-coordinate  defined  here  as  the  coordinate  along  the  waveguide  axis . 

The  preliminary  study  was  restricted  to  homogeneous  waveguides  whose  cross- 
section  and  material  properties  cure  z-independent . Furthermore,  as  already 
stated  the  analysis  was  limited  to  waveguides  whose  transverse  dimension 
was  small  enough  to  eliminate  thickness  vibrations,  e.g.  Euler  beams  on 
distributed  springs,  simply-supported  strips,  and  thin  cylindrical  shells. 

In  other  words,  three-dimensional  elasticity  theory  was  not  required.  In  I 
the  next  section,  we  shall  review  the  current  status  of  the  theory  of 
uniform  three-dimensional  solid  elastic  waveguides.  In  Section  3,  new 
theoretical  results  generated  from  the  present  study  will  be  stated  for 
three-dimensional  semi- infinite  z-independent  waveguides.  In  Section  4, 
the  properties  of  the  fundamental  mode  of  propagation  in  semi-infinite  wave- 
guides with  z-dependent  parameters  are  discussed.  Starting  with  a 
review  of  the  existing  literature,  we  proceed  to  analyze  non-uniform  wave- 
guides which  have  apparently  not  been  analyzed,  particularly  waveguides 
with  rapidly  varying  physical  constants  and  cross  sections. 

Results  generated  in  this  study  and  solutions  to  problems  which  have 
to  our  knowledge  not  been  published  are: 

1.  In  three-dimensional  elastic  cylinders  conducting  longi- 
tudinal axi-symmetric  waves,  the  axial  modal  impedance  defined  as  the  ratio 
of  the  axial  stress  averaged  over  the  cross  section,  divided  by  the  similarly 

averaged  social  velocity  equals  pc  , where  c is  the  modal  phase  velocity. 

z z 

The  simplicity  of  this  result,  reminiscent  of  acoustic  waveguides,  is 
surprising,  because  both  a dilatational  smd  a shear  potential  are  required 
to  describe  the  wave  motion. 
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2.  For  the  same  waveguide,  an  expression  was  obtained  for  the 
radial  modal  impedance,  defined  as  the  ratio  of  the  radial  shear  stress 
averaged  over  the  cross  section  divided  by  the  radial  velocity  of  the 
cylindrical  surface. 

3.  At  cut-off,  the  radial  modal  impedance,  and  consequently 
the  resultant  radial  drive-point  impedance  vanishes.  Simultaneously,  the 
axial  modal  impedance  diverges. 

4.  In  non-uniform  waveguides  located  in  the  region  z > 0,  whose 
properties  change  either  exponentially  or  as  a power  of  z,  the  effect,  on 
the  decay  of  the  fundamental  mode  with  increasing  z,  of  a change  in  cross 

section  has  precisely  twice  the  magnitude  of  a change  in  density  or  in  the 

* 

Young's  modulus  (see  Table  1). 

5.  The  phase  velocity  of  the  fundamental  mode  equals  the  local 
"bar"  velocity. 

6.  Exceptionally,  viz.  for  a constant  phase  velocity,  the 
fundamental  mode  of  non-uniform  waveguides  with  exponentially  varying 
parameters  displays  a cut-off  frequency.  This  is  a generalization  of  the 
well  known  property  of  acoustic  horns. 

For  the  non-uniform  waveguides  conducting  longitudinal  waves,  analytical 
solutions  have  been  constructed  in  the  long  wavelength  limit,  for  the 
following  situations  (see  Table  1) : 

a.  The  density.  Young's  modulus,  and  cross  sectional  area 
vary  exponentially  with  the  axial  coordinate  z at  different  rates,  no 
restriction  being  placed  on  the  rate  of  change  with  z except  in  so  far  as 
radial  displacement  are  not  accounted  for.  This  applies  to  subsequent 
configurations  as  weJl. 

b.  The  density.  Young's  modulus,  and  cross  sectional  area 
each  vary  as  a different  arbitrary  power  of  z. 

* This  conclusion  is  based  on  the  asymptotic  large-argument  or  far-field  form 
of  the  solutions  in  Eqs.  29  and  36.  It  was  pointed  out  to  the  authors  that 
this  asymptotic  expression  is  tantamount  to  the  WKB  approximation  which  does 
place  some  restriction  on  the  rate  of  variation  of  the  waveguide  parameters. 
Our  conclusions  regarding  the  respective  effects  of  E,  p,  and  S are  not 
restricted  to  the  two  analytically  tractable  classes  of  waveguides  analyzed 
here  but  can  be  generalized  to  other  waveguides  to  which  the  WKB  approximation 
is  applicable. 
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7.  For  flexural  waveguides  for  which  the  product  of  bar 
1/2  1/2 

velocity  (E/p)  times  radius  of  gyration  (I/S)  varies  parabolically  or 
linearly  with  z,  e.g.  a cone  or  a wedge,  the  long-wavelength  solution  for 
the  fundamental  antisymmetric  (i.e.  flexural)  mode  was  constructed,  both  for 
the  free  waveguide  and  the  waveguide  supported  on  a distributed  spring  whose 
stiffness  varies  like  the  cross  sectional  area.  The  spring-supported  wave- 
guide displays  a cut-off  frequency  at  the  natural  frequency  corresponding  to 
rigid-body  translational  vibration. 

B.  Dispersion  and  Cut-off  Characteristics  of  Three -d -i  m^nsional 
Elastic  Waveguides 

1.  Torsional  Waveguides 

The  theory  of  torsional  wave  motions  in  cylindrical  wave- 
guides, as  well  as  that  of  congressional  and  flexural  waveguides,  is 

2 3 4 

presented  in  a number  of  classical  texts.  ' ' It  need  therefore  not  be 
paraphrased  here.  The  torsional  waveguide  differs  from  the  latter  two  in 
that  it  requires  only  one  potential  to  describe  its  dynamic  behavior.  In 
contrast,  compressional  and  flexural  waveguides  require  both  a dilatational 
and  a shear  potential  as  will  be  illustrated  in  Section  C. 

The  solution  of  the  torsional  waveguide  is  as  straight-forward  as 
that  of  an  acoustic  waveguide,  the  single  potential  required  to  describe 
its  response  being  the  solution  of  the  wave  equation  or,  for  the  harmonic 
conditions  assumed,  of  the  Helmholtz  equation.  It  is  therefore  a simple 
matter  to  match  an  excitation,  viz.  an  oscillatory  torque,  to  the  modal 
series  describing  the  prescribed  shear  stress  distribution  over  the  cross 

5 

section  being  excited.  Onoe  was  thus  able  to  construct  the  driving  point 
impedance  or,  in  his  formulation,  its  reciprocal,  the  admittance,  as  a 
modal  series,  in  the  same  simple  form  arrived  at  for  the  elementary 
structural  waveguides  analyzed  in  Ref.  1.  As  in  the  latter  study,  Onoe’s 
drive-point  admittance,  defined  as  the  angular  velocity  divided  by  the 
oscillatory  torque,  becomes  infinite  at  every  cut-off  frequency.  As 
explained  in  Ref.  1,  the  reason  is  that  the  admittance  is  a modal  series 
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whose  terns  have  a denominator  proportional  to  the  axial  wavenumber  of 
the  various  inodes.  By  definition,  this  wavenumber  vanishes  at  the  cut- 
off frequency  of  a given  mode.  This  is  apparent  if  one  examines  Onoe’s 
Eq.  29,  and  notes  that  SL  is  an  axial  wavenumber.  The  reason  the  resultant 
impedance  vanishes  (i.e.  the  resultant  admittance  becomes  infinite)  is 
amply  discussed  in  Ref.  1 and  need  not  be  repeated  here. 

2.  Compressional  and  Flexural  Waveguides 

For  these  waveguides  dispersion  curves  have  been  plotted  and 
impedances  have  occasionally  been  computed  numerically.  This  may  be  an 
opportune  place  for  correcting  an  error  perpetuated  in  a familiar  mono- 
graph,  which  erroneously  states  that  the  characteristic  equation  of 
flexural  waves  in  cylindrical  waveguides  admits  only  one  root  corresponding 
to  a single  modal  dispersion  curve. 

The  reason  the  impedances  (or  admittances)  of  compressional  and  flexural 
waveguides,  have  not  been  formulated  rigorously  as  modal  series  is  that  the 
normal  modes  of  propagation  do  not  yield  orrthogonality  relations  suitable 
for  expressing  a prescribed  end  loading  as  a modal  series.5  Furthermore, 
the  Pochhammer-Chree  frequency  equation  which  governs  the  dispersion  curves 
of  cylindrical  waveguides  and  the  Rayleigh-Lamb  equation  which  governs  the 
dispersion  curves  of  plates,  admit  roots  in  the  form  of  complex  wave  numbers 
extending  down  to  vanishing  frequencies.  Other  roots  are  as  for  torsional 
or  acoustic  waveguides,  imaginary  below  cut-off  and  real  above.  The 
existence  of  complex  roots  was  first  pointed  out  by  Adem^  for  cylindrical 
waveguides  conducting  compressional  waves.  We  will  see  in  subsection  C.4 
that  the  values  of  drive-point  admittance  he  computed  from  a modal  series 
appear  to  be  incorrect  because  the  modes  do  not  form  a complete  orthogonal 
set.  The  existence  of  complex  wavenumbers  has  been  exhaustively  studied 
by  Onoe  et  al.  for  plates.  His  results  cure  reproduced  in  Ref.  4,  page  137. 
This  figure  shows,  in  particular,  the  existence  of  a large  number  of  higher- 
order  modes  which,  even  though  exponentially  damped,  display  a complex 
wavenumber,  i.e.  a finite  phase  velocity  at  vanishing  frequency,  incompatible 
with  the  cut-off  phenomenon.  Several  other  authors,  such  as  Mindlin  and 
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McNiven  have  obtained  similar'  results  both  for  cylinders  and  plates,  but 

there  is  no  need  here  for  aui  exhaustive  bibliography.  The  interested 

reader  is  referred  to  Miklowitz's  excellent  even  though  somewhat  dated 
8 

review  paper. 

Both  these  characteristics,  i.e.  the  unsuitability  of  the  normal 
modes  for  constructing  a series  expression  of  the  admittance,  and  the 
absence  of  true  cut-off  frequencies  for  a number  of  modes,  prevent  the 
analytical  formulation  of  the  drive-point  impedance  (see  subsection  C.4) . 

This  defeats  a portion  of  the  original  purpose  of  this  study.  We  earn, 
however,  formulate  some  interesting  results  concerning  modal  impedances 
not  formulated  by  earlier  authors  by  limiting  ourselves  to  inodes  displaying 
cut-off  frequencies,  i.e.  imaginary  or  real  rather  than  complex  wave 
numbers.  Even  though  the  drive-point  admittance  is  intractable  as  a modal 
series,  its  behavior  at  a modal  cut-off  frequency,  being  governed  by  a 
single  mode,  is  tractable. 

C.  Axial  and  Radial  Modal  Impedances  of  Longitudinal  Wave  Modes  in 
Solid  Cylinders 

1.  Statement  of  the  Characteristic  Equation,  Stresses,  and 
Displacements 

The  purpose  of  this  section  is  to  derive  some  novel  proper- 
ties of  the  modal  impedances  of  the  higher  modes  of  propagation  of  solid 
structural  waveguides  conducting  compressional  waves.  For  the  purpose 
of  this  analysis,  we  have  selected  the  axisymmetric  wave  motions  of  a 

semi-infinite  cylindrical  waveguide.  We  shall  start  from  the  equations 

3 9 

of  motion  as  formulated  by  Redwood,  page  137,  and  in  Davies'  monograph. 

For  longitudinal  excitation  the  results  will  be  expressed  in  the  form  of 

two  modal  impedances.  Hie  axial  modal  impedance  is  defined  as  the  ratio 

of  the  axial  stress  averaged  over  the  cross  section  divided  by  the  similarily 

averaged  axial  velocity.  The  radial  impedance  is  the  space-averaged 

axially  symmetric,  radially  oriented  shear  stress  divided  by  the  radial 

velocity  of  the  cylindrical  surface. 
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The  displacements  are  derived  from  a dilatational  potential  <j>  and  a 
shear  potential  ip.  For  axi symmetric  wave  motions,  both  potentials  are 
governed  by  the  two-dimensional  Helmholtz  equation  in  z and  r: 


v2<p  + ~ ♦ - 0 


v2^  + ~ ip  = 0 

ct 


(1) 


Here  and  are  respectively  the  velocity  of  dilatational  and  transverse 

waves  in  the  elastic  medium.  If  X and  y are  Lame's  constants  and  p the 

1/2  1/2 

density,  the  former  velocity  equals  [(X  + 2y)/p]  and  the  latter  (y/p) 
Selecting  solutions  of  the  form 

<b  - <J>Q(r)  expfiyz-iwt) 


and  similarly  for  ipQ,  one  finds  that  and  ipQ  are  Bessel  functions 

2 to2  2 

(r)  = A J (hr)  expdyz-iut) , h = — r - y 

O fc 

cd 

(2) 

2 a)2  2 

ip(r)  - C J (kr)  exp(iyz-io)t) , k * -j  - y 

ct 

Note  that  <P  and  A have  units  of  length  squared,  while  ip  and  C have  units  of 
length  cubed.  The  corresponding  radial  and  axial  displacements  are 

u = 1±  + 1jL_ 

r dr  9r9z 


_ 9_£  _ 9jj>_  _ 1_ 

Uz  9 z .2  r 9r 
9r 


(3) 
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This  impedance  governs  the  axial  modal  response  to  an  axial  excitation. 
It  is  defined  here  as 


Z 

z 


a 

{ t (r)  rdr 
0 zz 

a 

/ u (r)  rdr 
0 2 


The  integrals  are  readily  evaluated  by  noting  that 


(8) 


a 

f 

0 


J (hr)  rdr  = ~ J,  (ha) 
o hi 


(9) 


and  similarly  for  J (kr) . Consequently,  the  integral  in  the  numerator 
yields 


a 


T rdi 
zz 


. 2 2 
+ hj*. 
2 .2 
Y -k 


Jj_  (ha) 


(10) 


2 

Multiplying  through  by  cfc,  reducing  the  terms  in  brackets  to  the  same 
denominator,  and  noting  that,  from  Eq.  2, 


J S 


if 


I 

i 

I 


2 <ji 

°t  .2  2 

k +y 

Eq.  10  becomes 

2 

a Apoj  aJ,  (ha) 

- / t rdr  = — 

0 22  h(Y~-k~) 

The  denominator  of  Eq.  8 similarly  yields 


j ma;  ( 2 2 2\ 

-J-  (T  ■*  «*  ) 


* . . Atoya  ( , . 2h^ 

{uzrdr=  h (1  + 7t:2- 

0 \ (y  -k  ) 


Jx(ha) 


Amya  J (ha)  _ ^ 7 

r— (y  ”k“+2n  ) 

h(Y-k) 
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Substituting  Eqs.  11  and  12  in  Eq.  8,  the  impedance  now  takes  the  simple 
form 

Zz  = pw/Y  (13) 

Referring  to  Eq.  2,  we  note  that  the  axial  wavenumber  is  related  to  the 
phase  velocity  as 

Y = oj/c 

z 

The  impedance  thus  reduces  to 

Z * pc  (14) 

z z 

The  axial  modal  impedance  of  the  compress  ional  waveguide  is  in  the  form  of 
the  plane-wave  impedance . This  result  is  remarkable  by  its  simplicity  if 
one  considers  that  two  potentials  characterized  by  two  different  charac- 
teristic velocities  c and  c are  required  to  formulate  the  modal  impedance . 
9 a 

Even  though  Davies  evaluates  this  impedance  his  eaqpression  is  lengthy  and 

involved,  being  expressed  in  terms  of  auxiliary  functions.  He  therefore 

did  not  realize  the  potential  compactness  and  the  physical  meaning  of  this 

impedance.  Since  cut-off  is  characterized  by  y=0,  c =°°,  Eqs.  13  and  14 

z 

indicate  that  the  axial  modal  impedance  becomes  infinite  at  cut-off. 


3.  The  Radial  Modal  Impedance 

This  impedance  is  defined  here  as  the  radially-oriented 
shear  stress  averaged  over  the  cross  sectional  area  divided  by  the  radial 
velocity  u^fa)  on  rhe  cylinder  surface. 


Z = 
r 


2f  x rdr 

o zr 

2.  , , 

a u^(a) 


(15) 


The  radial  velocity  is  readily  constructed  from  Eqs.  2 and  3, 


u (r)  = imhA 
r 


2 2 Jx(ha) 

2,2  J,  (ka)  Ji(kr)  " ,Tl(hr) 


LY  -k 


(16) 


i 


L 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC. 


Setting  r=a,  this  reduces  to 


. . . i»hMY2+k2)  ....  “ “ Jl(ha) 

Ur  2j~  • Jl(ha>  * 2,  2 2. 

Y -*  c^(y  -k  ) 


(17) 


where  use  has  been  made  of  Eq.  10a.  The  shear  stress  is,  from  Eq.  4c, 


zr 


2 T Jl(ha) 

= -i2PctYhA  [j^hr)  - Jx 


(kr) 


(13) 


. 10 


The  integrals  in  Eq.  15  can  be  evaluated  by  means  of  the  relation 
a 


f J..  (hr)  rdr  ~ |£  [J. (ha)  (ha)  - H (ha)J  (ha)] 
0 l zn  l ~o  1 o 


Where  is  the  Struve  function.  Substituting  this,  and  the  equivalent 


egression  for  the  integral  of  J^kr)  , the  numerator  in  Eq.  15  becomes 


21  T^rdr  = - i2irpctAay  ^(ha) 


, 1U  (ka)H-(ka) 

H0(ha)  --H0tta)  + 


- J (ha)H. (ha) 
o 


(19) 


The  radial  impedance  finally  becomes 


Z = 
r 


4 2 2 

2irpctY(Y  -k  ) 


r 


ahuf 


-H0(ha)  +- 


HQ(ka)  - 


J (ka)  H (ka) 
o ~1 


J-^ka) 


J (ha)H,  (ha) 
o ~1 


(ha, 


(20) 


/ 

At  cut-off,  the  modal  radial  impedance  vanishes,  since  y=0*  The  cut-off 
frequency  is  therefore  seen  once  again  to  represent  a resonance  of  the  two- 
dimensional  z- independent  standing  wave  system  of  the  waveguide  cross 
section.  It  was  shown  in  Ref.  1 that  the  divergence  of  a single  mode  is 
sufficient  to  cause  the  drive-point  admittance  to  diverge. 
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4.  Discussion  of  the  Drive-point  Impedance 

If  the  stresses  and  displacements  given  above  were  functions 

of  a single  potential,  the  boundary  condition  prescribed  over  the  driven 

cross  section,  e.g.  t (r)  could  be  simply  computed  by  expanding  it  as  a 

zz 

series  of  Bessel  functions,  and  using  the  orthogonality  relation  of  Bessel 
functions  to  compute  the  modal  coefficients.  More  generally,  this  method 
can  be  used  when  the  boundary  condition  takes  the  form 

J (ka)  + Aj’(ka)  = 0 (21) 

n n 

where  A is  a constant  and  k a wavenumber.  The  fact  that  t , Eq.  7 contains 

zz 

two  Bessel  functions,  one  for  each  potential,  eliminates  this  approach  even 
though  each  potential  separately,  being  a solution  of  the  Sturm-Liouville  equa- 
tion can  be  expanded  in  an  orthogonal  set  of  eigen  functions . ^ Because  their 
stresses  are  expressible  with  a single  potential  satisfying  that  equation, 
torsional  waveguides  are  the  one  3-dimensional  elastic  waveguide  whose  impedance 
is  analytically  tractable.5  The  coupling  of  eigen  modes  in  thick  cylindrical 

shells  rigorously  analyzed  as  a problem  in  three-dimensional  elasticity  is 

12 

noted  by  Armenakas  et  al.  We  must  be  uneasy  about  the  manner  in  which 
7 

Adem  formulates  the  forced  response  of  a cylindrical  waveguide  m terms  of 
its  eigen  functions:  "If  we  know  only  some  of  the  roots  (of  the  character- 
istic equation),  then  for  the  other  roots  we  can  take  B(£^)  =0,"  where 
B(£^)  is  the  amplitude  of  the  eigen  mode  of  order  q.  Quite  clearly  the 
fact  that  one  can  arbitrarily  eliminate  some  eigen  functions,  shows  that 
they  do  not  form  a complete  orthogonal  set. 

The  orthogonality  of  the  modes  in  infinite  or  semi-infinite  places 

13 

(Lamb  waves)  was  examined  by  Lyon.  He  was  apparently  the  first  worker 
to  point  out  that  these  modes  do  not  form  an  orthogonal  set  with  respect 
to  the  plate  thickness.  This  problem  has  been  receiving  some  attention  in 
the  Russian  literature.  ' Even  though  the  latter  two  papers  construct 

some  orthogonality  relations,  they  are  not  suitable  for  the  series 
expansion  of  applied  loads,  i.e.  for  the  construction  of  the  drive-point 
impedance  as  a modal  series. 
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D.  The  Fundamental  Mode  of  Non-uniform  Waveguides  Conducting 
Longitudinal  and  Transverse  Waves 

1.  General  Discussion 

The  motions  in  axi symmetric  waveguides  of  exponentially 
varying  cross  section,  i.e.  in  solid  horns,  have  been  studied  by  means  of 
an  approximate  theory  valid  when  the  transverse  dimension  is  small  in 
terms  of  wavelengths,1^  thus  automatically  eliminating  three-dimensional 
cut-off  phenomena. 

Xeller  and  one  of  his  students  constructed  an  asymptotic  WKB-type 

solution  for  higher  modes  including  flexural,  torsional,  and  longitudinal 

wave  motions,  restricted  to  solid  waveguides  with  a slowly  varying  cross 
17 

section.  This  restriction  results  in  a solution  which  predicts  local 
z-dependent  cut-off  frequencies  equal  to  those  of  the  uniform  cylinder 
whose  local  cross  section  equals  that  of  the  non-uniform  waveguide  cross 
section  at  z.  In  view  of  the  difficulties  mentioned  in  the  preceeding 
section,  it  is  not  surprising  that  the  authors  made  no  attempt  to  evaluate 
drive-point  impedances. 

We  can  however,  draw  certain  conclusions  from  Ref.  17,  We  have  seen 

that  the  axial  modal  impedance  becomes  infinite  at  cut-off.  Consequently, 

if  we  consider  acoustic  energy  propagating  in  the  direction  of  decreasing 

cross  sections,  all  modes  but  the  plane-wave  mode  will  gradually  be 

reflected  back,  thus  giving  rise  to  standing  waves.  A similar  conclusion 

18 

is  reached  for  acoustic  waveguides  of  slowly  varying  cross  section. 

We  shall  examine  some  specific  non-uniform  waveguides  which  admit  an 
analytical  solution  but  which  have  not  yet  been  treated  in  the  literature. 
We  shall  obtain  dispersion  relations  for  the  fundamental  modes  of  longi- 
tudinal and  flexural  waveguides,  but  shall  not  construct  the  Green's 
influence  functions. 
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2.  Longitudinal  Waves  in  Waveguides  with  Density,  Young’s  Modulus, 


and  Cross-section  Varying  Exponentially  at  Arbitra 


Rates 


These  quantities  vary  as 


E = E e 
o 


S * S e 
o 

The  constants  a,  e,  and  6 are  not  restricted  in  magnitude.  For  the  long- 
wavelength  situation  exclusively  considered  throughout  this  Section  D,  radial 
motion  associated  with  the  fundamental  mode  can  be  neglected.  The  steady- 
state  equation  of  motion  of  a non-uniform  waveguide  is  therefore  an  ordinary 
D.E.  formulated  in  terms  of  the  axial  displacement,  which  need  no  longer  be 
identified  by  the  subscript  z: 


d du.  , 2 „ _ 

— (ES  — ) + «o  pSu  - 0 


Differentiating,  and  dividing  through  by  E S,  one  obtains  the  equation 


d2u  / .1  dE  + 1 dS  \ du 

2 \E  dz  S dz ) dz 


+ “_£  u = o 
E 


For  the  parameters  described  in  Eq.  22 
2 

+ (e-*6  ) + k2  exp  [ (o-e)  z]  u = 0 

dz 


where  we  have  used  the  bar  velocity  and  wavenumber 


cq  = (E0/P0)  » c(z)  = cq  exp[(£-a)z/2] 


kQ  = tn/cQ  , y(z)  = kQ  exp[(a-e)  z/2] 
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I 


1 

5 I 


We  now  make  the  transformation  of  variables 

u = u exp( — z) 

du  / du  e-:d  -A  ( e-*d  \ 

dz  ~ \dz  ~ 2 / ^ 2 / 

d2u  d2u  , ,,  du  . (eH6  - / e-»6  \ 

|v"  5S  + -V- UJ  -* 


Substituting  in  Eq.  25/  one  obtains  the  equation  governing  u: 


>- 

i* 

s L 


.2  , (e+S) 

k exp(a-e) z — 

O 4 


u = 0 


We  now  perform  a transformation  in  the  independent  variable 

2 - 
Z = Z 

ot-e 

The  equation  now  becomes 


2-  f 4k2 
i u + o 

2 . (ct+e) 


2z  (e-6 ) ' 

a — — ; 

(o-e) ' 


This  in  the  form  of  Bessel's  equation.  The  transformed  variable  u is 
obtained  in  the  form  of  a cylinder  function  C of  order  (e-hS ) /(a-e)  and 
argument  2kQe^/(a-e) . Transforming  back  to  the  original  variables,  we 
finally  have 


u(z)  = U exp 


The  particular  cylinder  function  selected  must  be  well  behaved  in  the  region 
of  z relevant  to  the  situation  being  investigated.  We  shall  construct 
general  solutions  but  shall  not  illustrate  the  matching  of  linear  combinations 


||  of  cylinder  functions  to  particular  excitations.  We  shall  merely  concern 

it  ’ ourselves  with  the  dispersion  characteristics  of  the  phase  velocity.  For  this 

fr 

jj|.  * purpose,  we  consider  waves  propagating  in  the  positive  z-direction.  Assuming, 
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e+S  >0,  a-e>0,a  suitable  solution,  i.e.  one  which  converges  as 
z+*>,  is  the  Hankel  function  of  the  first  kind..  For  large  argument 

r - 


H 


(1) 


1/2 


/a-e  \ . i^o 

(«*>/<«-.>  = liiTj  “W  z)  sxpi  ' 


(30) 


exp 


where  v is  the  order  of  the  Hankel  function.  When  this  is  substituted  in 
Eq.  29,  the  solution  becomes 

1/2 

. ) 

o 


, . , „ / a-e  \ ' l a+e+25  \ id> 

uU,t>  - a j exp  (-  — j—  z)  e 


(31) 


where 


2k 

o 

<t>  = — exp 
a-e 


03t  - £ (2v+l) 
4 


(32) 


The  variation  of  the  amplitude  of  the  solution  with  increasing  z is  embodied 
in  the  real  exponential  term.  The  change  in  cross  sectional  has  twice  the 
effect  of  either  the  density  or  the  Young's  modulus  for  comparable  exponential 
coefficients.  Signal  attenuation  by  the  increase  in  cross  section  is  in  the 
nature  of  a spreading  loss.  The  attenuation  associated  with  the  increase 
in  Young's  modulus  and  density  embodies  gradual  backscattering  as  the  signal 
penetrates  into  a region  of  ever  increasing  characteristic  impedance.  For 
negative  powers  of  the  exponentials  in  Eq.  22,  the  signal  increases,  as 
anticipated.  The  solution  displays  z-dependence  on  the  coefficient  deter- 
mining the  cross  sectional  area  variation  that  is  precisely  twice  that  of 
either  modulus  of  density  coefficients.  The  phase  velocity  is  obtained  from 
the  phase  angle 


— 3 <j>  9 z 

C “ " 9t  9<j> 


(33) 


where 


dd> 

9t  = 


■w 


94>  , 

r-1-  = k exp 
9z  o c 


(¥■) 


(33) 
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We  can  solve  for  a dispersive  phase  velocity  which  equals  the  local  bar 
velocity,  Eq.  26. 


When  a-*e,  i.e.  when  the  bar  velocity  is  z- independent,  the  order  of 
the  Hankel  function  diverges,  as  does  its  argument.  Returning  to  the 
original  differential  , Eq.  25,  we  note  that  the  coefficient  of 

the  linear  term  becomes  z-independent . The  solution  is  an  exponential, 
say  ejq>(Nz) , where  N is  a solution  of  the  characteristic  equation 

N2  + (e+6)N  + k2  = 0 
o 


This  is  a quadratic  equation  with  two  roots.  The  axial  displacement  now 
becomes 


u(z)  = 0 exp 


a = e (34) 


This  is  an  exponentially  damped  solution  which  admits 
waves  in  the  frequency  range  where  the  square  root  is 
the  cut-off  frequency. 

(e-*6  )c 


co 


4ir 


a = e 


damped  propagating 
imaginary,  i.e.  above 


(35) 


If  a = e = C.  this  cut-off  frequency  reduces,  as  anticipated,  to  the 
familiar  result  obtained  for  the  exponential  horn. 


3.  Longitudinal  Waveguide  Whose  Density,  Young's  Modulus,  and 
Cross  Sectional  Area  Vary  as  Arbitrary  Powers  of  z 

The  density.  Young's  modulus,  and  cross  sectional  area, 

vary  as 


* n 
p = Az 


= - <B/A)1/2  z<™)/2 
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Eg.  24  becomes 


2 2 
d u . m+r  du  , m A n-m  n 

— — + — + — — - z u = 0 

,2  z dz  B 
dz 


20 


This  is  a form  of  Bessel's  equation, 
cylinder  function  of  fractional  order: 


. . _ (1-m-r) /2 

u(z)  = Uz  C 


(35) 


Its  solution  is  proportional  to  a 


(1-m-r) /(n-m+2) 


2uAV22(n-«2)/2 


(n-m+2) B 


1/2 


(36) 


As  in  the  discussion  of  Eg.  29,  a Hankel  function  of  the  first  kind  is 

selected  to  represent  waves  propagating  in  the  positive  z-direction.  The 

resulting  phase  angle,  when  substituted  in  Eq.  33,  leads  to  the  usual 

conclusion  that  the  phase  velocity  equals  the  local  bar  velocity,  Eq.  34. 

In  this  same  region,  the  cylinder  function  converges  as  z m+2^4.  The 

amplitude  of  the  plane  wave  therefore  converges  as  z (n+m+2r)/4>  we  note 

that  as  for  exponentially  varying  waveguide  parameters,  Eq.  31  an  increase 

in  density  and  in  Young's  modulus  bring  about  a conparable  acceleration  in 

convergence  of  the  solution  with  increasing  z,  and  that  a change  in  cross 

section  has  an  effect  of  precisely  twice  this  magnitude.  Negative  powers 

of  z in  Eq.  34  produce,  as  anticipated,  an  enhancement  of  the  signal.  The 

solution  in  Eq.  36  ;'s  expressible  in  terms  of  familiar  cylinder  fractions 

for  specified  z-dependences  of  the  waveguide  parameters.  The  results  are 

listed  on  Table  1.  The  cases  which  admit  a solution  in  form  of  Airy 

fractions  are  related  to  similar  solutions  for  acoustic  waveguides  with 

18 

variable  cross  sections  and  for  acoustic  waveguides  with  sound  velocity 

21 

gradients  and  absorptive  boundaries . Both  these  analyses  of  acoustic 
waveguides  were  carried  through  in  sufficient  depth  to  include  higher  modes 
and  cut-off  frequencies.  The  greater  complexity  of  the  solid  elastic  wave- 
guide requires  a laborious  analysis  probably  not  justified  by  the  limited 
practical  interest  of  this  waveguide  configuration.  When  m=n=0,  the  solution 
reduces  to  the  familiar  solution  for  acoustical  horns.  The  most  widely  used 
is  the  conical  horn  (r=2) , which  does  not  display  a cut-off  frequency  in 
contrast  to  the  exponential  horn  mentioned  in  the  preceding  section. 


WS£WBWCT1B5W ' ■raWL_»BE; 
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When  (m-.  )=>2,  Eq.  36  becomes  indeterminate.  The  differential  equation 
of  motion  becomes 


2 2 
d u m+r  du  to  A _ 

, 2 + z dz„2U“ 
dz  Bz 


m-n=2 


M 


The  solution  is  in  the  form  of  U/z  . Substituting  in  Eq.  37,  one  constructs 
the  characteristic  equation 

2 2 

M -M(nH-r-l)  + to  (A/B)  = 0 , m-n=2 

Solving  this  quadratic  equation,  one  obtains  a solution  in  the  form 

(m+r-1)  - f (m+r-1) 2 


u(z)  = exp 


(-  | iatEii  ; ^i£±£zlli  . iL* ] 1/2  | 2)  ’ 


m-n  = 2 


(38) 


This  deformation  is  in  the  form  of  a non-propagating  near-field  at  small 
frequencies  where  the  exponent  is  real.  The  response  is  a damped  propagating 


wave  when  the  power  is  coup  lex,  i.e.,  in  the  frequency  range 

1 / B \ 1/2  m+r-1 

2 \a)  2 


f > 


To  conclude  this  section,  we  turn  our  attention  to  the  fundamental 
mode  of  non-uniform  flexural  waveguides. 


4.  Non-uniform  Flexural  Waveguides 

In  the  low-frequency  limit,  the  Bernoulli-Euler  equation 
which  governs  the  steady-state  flexural  vibrations  of  non-uniform  beams  is 


.4  A2  a2  o 

El  — ~ + — - (El)  — - - pSu  w = 0 
dz  dz  dz 


This  takes  the  more  explicit  form 


(39) 


d4w  / 1_  d2I  1 d2E  _2_  dE  dl  \ d2w  pSaj2w  _ 

J4+\I  ,2+E  ,2  EIdzdz  ,2'  El 

dz  \ dz  dz  / dz 


(40) 


For  a uniform  beam,  the  coefficient  of  the  second  derivative  vanishes,  and 
that  of  the  linear  term  is  constant.  The  solution  takes  on  the  familiar 
form 
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w(z)  = W e^2 


where  y is  the  flexural  wavenumber 


Y * (k/r) 
o g 


where 


cq  = (E/p) 


r - (I/S)-1-^ 

9 

2 2 

Eq.  40  reduces  to  Bessel's  equation  if  the  d w/dz  -term  vanishes  and  if 

_2 

the  coefficient  of  the  linear  term  is  proportional  to  z . Both  these 

-2 

conditions  are  satisfied  if  the  product  pS  varies  as  z , E and  I being 
constant.  This  is  of  course  an  unrealistic  assumption.  We  can  however 
construct  a meaningful  mathematical  model  by  taking  E,  p,  and  hence  c^ 
constant,  and 


r = Bz 

g 


Y - (kQ/Bz) 


This  condition  is  satisfied  by  a wedge , for  which 

S = Az,  I = AB2z3  ( 

and  a cone  or  pyramid,  for  which 
2 2 4 

S = Az  , I = AB  z i 

Furthermore,  to  make  the  D.E.  tractable,  the  coefficient  of  the  second 


(43a) 


(43b) 


derivative  in  Eq.  40  must  be  negligible.  The  solution  thus  obtained  will 
be  shown  to  be  in  the  form  of  a cylinder  function.  In  the  large  argument 
limit,  this  solution  yields  the  second  derivative 


d2w  2 

— - = -4y  w 

dz 


yz  » 1 


The  second  derivatives  cam  be  ignored  if  y z « 24  for  the  wedge,  or  « 48 
for  the  cone,  conditions  obviously  not  satisfied  in  the  region  adjoining  the 
waveguide  apex,  where  z is  small  in  terms  of  wavelengths. 
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The  approximate  differential  equation  is  now  formally  similar  to 

the  flexural  equation  of  motion  of  uniform  beams,  even  though  the  coefficient 

-2 

of  the  linear  term  varies  as  z : 


d4w  4 

—4  “ Y w = u 

dz 


«■>  k 
2 o 


(44) 


This  equation  admits  a solution  in  the  form  of  a cylinder  function 


22 


w = WzC2(2iqyz)  = WzC2  [2iq(kQz/B) *^2] , q = 0,  1,  2 or  3 


(45) 


It  is  interesting  to  note  that  Kirchhoff  solved  the  problem  of  the  finite 
wedge  in  terms  of  cylinder  functions  as  early  as  1879. 

Eq.  45  embodies  propagating  waves  in  the  form  of  Hankel  functions  of 
order  2.  For  large  argument,  the  function  is  proportional  to  exp(i<j>)  where 

<j>  - 2(kQz/B)1/2-u)t 


The  corresponding  phase  velocity  computed  from  Eq.  33  is 

c = (wBz)  ^2  = (oic  r ) ^ 2 
o g 

As  usual,  the  phase  velocity  at  z is  the  phase  velocity  which  would  be 
observed  in  a uniform  waveguide  whose  material  properties  and  cross  section 
coincide  with  those  of  cross  section  z. 

To  conclude  this  section  we  consider  a non-uniform  flexural  waveguide 
whose  fundamental  mode  displays  cut-off  behavior.  For  a waveguide  mounted 
on  a distributed  spring,  the  equation  of  motion  becomes 

d4  2 

(El  — 7 + K - (1)  pS)  w = 0 
dz 


or 


2 

0) 


2 2 


c r 
o g 


w = 0 


(46) 
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where  K is  the  spring  stiffness  per  unit  length.  The  waveguide's  z-dependence 
is  still  governed  by  Eqs.  43.  The  spring  stiffness  is  selected  to  vary 
linearly  with  z for  a wedge.  Eg.  43a, 

K = Dz 

and  parabolically  for  a cone  or  pyramid,  Eq.  43b, 

K = Dz2 


For  either  configuration  the  ratio  K/pS  is  z- independent.  It  is  in  fact  the 
natural  frequency  o>n  of  the  undeformed  beam  undergoing  translational,  i.e. 
z-independent  vibrations.  Using  Eq.  44  we  have 


The  solution  of  this  equation  is  of  the  form  of  Eqs.  43  and  44,  but  the 

2 2 1/4 

flexural  wave  number  is  multiplied  by  (1  - a>  /o>  ) 

n 

The  construction  of  the  dispersion  curves  proceeds  as  in  the  preceding 
subsection: 


With  increasing  frequency,  this  phase  velocity  tends  to  Eq.  45.  It  diverges 
as  and  is  complex  in  the  frequency  range  co <^n-  The  wave  therefore 

attenuates  exponentially  as  it  propagates.  The  frequency  to  = u>n  is 
associated  with  a zero  drive-point  impedance,  a transverse  force  exciting 
the  rigid-body  resonance.  Like  the  phase  velocity,  the  characteristic 
impedance  of  flexural  waves  is  infinite  at  this  frequency. 

This  cut-off  phenomenon  is  observed  in  uniform  spring-mounted  flexural 
waveguides.^-  The  condition  for  its  occurrence  in  a waveguide  with  variable 
mass  per  unit  length  is  that  the  spring  stiffness  and  the  mass  of  the  beam 
display  the  same  z-dependence,  thus  permitting  rigid-body  translational 
vibrations  uncoupled  from  any  rotational  motion. 


i 
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Table  1 - Effect  of  Various  Parameters  on  the  Propagation  of  the 
Fundamental  Mode  in  Non-uniform  "Longitudinal"  Waveguides 
(All  Coefficients  are  Positive;  the  Waveguide  Extends  over  the  Region  z>0) 
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